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Abstract. We investigate the hydrodynamics of three- 
dimensional classical Bondi-Hoyle accretion. A totally absorb- 
ing sphere of different sizes (1, 0.1 and 0.02 accretion radii) 
moves at different Mach numbers (0.6, 1.4, 3.0 and 10) rela- 
tive to a homogeneous and slightly perturbed medium, which 
is taken to be an ideal, nearly isothermal, gas (7 = 1.01). We 
examine the influence of Mach number of the flow and size of 



the accretor upon the physical behaviour of the flow and the ac- 
cretion rates. The hydrodynamics is modeled by the "Piecewise 
^ — I Parabolic Method" (PPM). The resolution in the vicinity of the 
accretor is increased by multiply nesting several 32 3 -zone grids 



around the sphere, each finer grid being a factor of two smaller 



v "-^ in zone dimension than the next coarser grid. This allows us to 
include a coarse model for the surface of the accretor (vacuum 
sphere) on the finest grid while at the same time evolving the 
t-H gas on the coarser grids. 

For small Mach numbers (0.6 and 1.4) the flow patterns 
^ tend towards a steady state, while in the case of supersonic flow 
(Mach 3 and 10) and small enough accretors, (radius of 0.1 
and 0.02 accretion radii) an unstable Mach cone develops, de- 
Q\ stroying axisymmetry. The shock cones in the supersonic mod- 
els never clear the surface of the accretors (they are tail shocks, 
O pot bow shocks) and the opening angle is smaller (compared 
' to models with larger 7) especially for the highly supersonic 
5^ models. The densities in the shock cone is larger for models 
^ with smaller 7. The fluctuations of the accretion rates and flow 
03 structures are weaker than in the corresponding models with 
larger 7. The hydrodynamic drag of all models with accretor 
sizes of 0.1 Ra or smaller acts in an accelerating direction, 
while the gravitational drag is always decelerating and larger 
than the hydrodynamic drag (thus the net force is decelerating). 
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1. Introduction 

This last instalment extends our investigation (started with Ruf- 
fert, 1994a, Ruffert & Arnett, 1994, Ruffert 1994b and 1995; 
henceforth referred to as Papers I to IV, respectively) of the clas- 
sic Bondi-Hoyle-Lyttleton accretion model. In this scenario a 
totally absorbing sphere moves with velocity relative to a 
surrounding homogeneous medium of density and sound 
speed Coo. One would like to find the accretion rates of vari- 
ous quantities (mass, angular momentum, etc) as well as the 
properties of the flow (e.g. distribution of matter and velocity, 
stability, etc). In this fifth instalment we investigate the accre- 
tion flow and accreted quantities for a nearly isothermal gas, 
which we modeled by taking the adiabatic index 7 to be 1.01 
and otherwise using the same parameters as had been done in 
the simulations of previous papers (which had 7 = 5/3 and 
7 = 4/3). In an astrophysical context, a gas might be accreted 
isothermally when the optical depth is small and the surplus 
energy (higher temperatures are generated due to gravitative 
compression of the medium) is radiated away. One of the sit- 
uations where this might be applicable is the accretion flow 
behind young stars moving relative to the ambient medium 
(M. Casali, private communication). Of course a simulation in- 
cluding radiative transport would be better, but it is numerically 
prohibitive. 

The Bondi-Hoyle accretion flow for specific heats 7 close 
to unity has been modeled previously by other authors, how- 
ever most of the simulations were two-dimensional. Shima et 
al. (1985) find that for 7=1.1 the shock is attached to the ac- 
cretor, the opening angle is smaller and the maximum densities 
are higher, than in cases with larger 7, and an accretion col- 
umn is formed. All these two-dimensional results have been re- 
produced and confirmed by our three-dimensional simulations. 
Petrich et al. (1989) extended the two-dimensional simulations 
also including a fully relativistic treatment. Within a factor of 
two, their results for the mass and momentum accretion rates 
coincide with the values of Shima et al. (1985) and from the few 
contour plots shown, also the density distribution seems similar. 
Matsuda et al. ( 1 992) performed several two-dimensional calcu- 
lations with 7 close to one and found that at first radial instability 



oscillations appear. The numerical resolution of the accretor in 
these simulations was, however, very crude. It was improved 
only slightly by Ishii et al. ( 1 993), who nonetheless attempted to 
additionally simulate three-dimensional isothermal flows. The 
latter authors conclude that also isothermal flows exhibit the 
flip-flop instability, but that in 3D it is much less violent than in 
2D. 

To more clearly delineate these differences, to compare the 
accretion rates of several quantities when flow parameters are 
changed and to obtain a more uniform picture of the accretion 
rates is the aim of this paper. In section 2 we give only a short 
summary of the numerical procedure used, since it has been 
described in the previous Papers I, II and III. Section 3 presents 
the results, which we analyze and interpret in Sect. 4. Section 5 
summarizes the implications of this work. 



2. Numerical Procedure and Initial Conditions 

Since the numerical procedures and initial conditions are mostly 
identical to what has already been described and used in Papers I 
to IV we will refrain from repeating every detail. Instead we 
only give a brief summary. 



2.1. Numerical Procedure 

A gravitating, totally absorbing "sphere" moves through an 
initially slightly perturbed (3%) but otherwise homogeneous 
medium. The relative velocity t>oo is varied in different models: 
we do simulations with Mach numbers of 0.6, 1.4, 3.0 
and 10. In the reference frame of the accretor the surrounding 
matter flows in +x-direction. Our units are (1) the sound speed 
Coo as velocity unit; (2) the accretion radius (Eq. 3) as unit of 
length, and (3) poo as density unit. Thus the unit of time is 
Ra/coc. 

The distribution of matter is discretised on multiply nested 
(e.g. Berger & Colella, 1989) equidistant Cartesian grids with 
zone size 6 and is evolved using the "Piecewise Parabolic 
Method" (PPM) of Colella & Woodward (1984). The equa- 
tion of state is that of a perfect gas with a specific heat ratio 
of 7 = 1.01. The model of the maximally accreting, vacuum 
sphere in a softened gravitational potential is summarized in 
Paper III. 



2.2. Models 

The combination of parameters that we varied, together with 
some results are summarised in Table 1 . The first letter in the 
model designation indicates the Mach number: E, F, G and H 
stand for Mach numbers of 0.6, 1.4, 3.0 and 10, respectively. 
The second letter specifies the size of the accretor: L (large), 
M (medium) and S (small) stand for accretor radii of 1, 0.1, 
and 0.02 R A , respectively. It is informative to compare the sizes 



analytically for spherically symmetric accretion (Bondi, 1952): 



d s = — — —Rb , 



with the Bondi radius given by 



Rb 



GM 



(1) 



(2) 



For further reference we define the accretion radius (Hoyle & 
Lyttleton, 1939, 1940a, 1940b, 1940c; Bondi & Hoyle, 1944) 



as 



Ra = 



2GM 



(3) 



in which M is the mass of the accretor, G the gravitational 
constant, the sound speed of the medium at infinity and 
the relative bulk flow of the medium at infinity. For a value of 
7 of 4/3 (as in Paper IV) the sonic distance is d s (7 = 4/3) = 
0.25 R B , while for 7 = 1.01 follows d s (j = 1.01) = 0.49i? B - 
Neglecting effects due to the relative motion of medium and 
accretor and assuming a roughly radial flow close to the accretor, 
one can boldly eliminate R B in Eq. 1 in favor of R A from Eq. 3 
and one obtains d s = R^M^iS — 3j)/i, where the Mach 
number is = i>oo/coo- Inserting both values of gamma, 

and all three Mach values of the supersonic flow (1.6, 3, 10) 
one finds, that the sonic distance is in all cases larger than the 
size of the medium (0. 1 R A ) and small accretor (0.02_R A )- Note, 
that this is only a rough estimate. 

The size of the largest grid L is adapted to the size of the 
accretor, larger accretors needing larger grids because the evo- 
lution time scales are longer and thus the gravitative influence 
reaches further out. The grids are nested to such a depth g that 
the radius of the accretor R spans several zones on the finest 
grid and the softening parameter e is then chosen to be a few 
zones less than the number of zones that the accretor spans. 

As far as computer resources permitted, we aimed at evolv- 
ing the models for at least as long as it takes the flow to move 
from the boundary to the position of the accretor which is at the 
center (crossing time scale). This time is given by L/2A4oo and 
ranges from about 1 to about 100 time units. The actual time tf 
that the model is run can be found in Table 1 . 

The calculations are performed on a Cray-YMP 4/64 and a 
Cray J90 8/512. They need about 12 MWords of main memory 
and take approximately 40 CPU-hours per simulated time unit 
(for the S = 1 /64 models and Mach 10; the 5 = 1 /256 models 
take four times as long, etc.; 6 is the size of a zone on the finest 
grid, see Table 1). 

3. Description of Results 

We will not only describe the results obtained in the new sim- 
ulations, but also compare these results with those obtained in 
Papers II, III and IV for models with 7 = 5/3 and 7 = 4/3. 
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Fig. 1. Contour plots showing snapshots of the density together with the flow pattern in a plane containing the center of the accretor for all 
subsonic models EL (panel a), EM (panel c), ES (panel e) and mildly supersonic models FL (panel b), FM (panel d), FS (panel f). The contour 
lines are spaced logarithmically in intervals of 0.05 dex for model EL and 0.1 dex for all other models. The bold contour levels correspond to 
log p = 0.01 for model EL and FL, log p = 0.2 for model EM, and log p = 1 .0 for models FM, ES and FS. The time of the snapshot together 
with the velocity scale is given in the legend in the upper right hand corner of each three panels. 
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Fig. 2. The accretion rates of several quantities are plotted as a function of time for all subsonic models EL, EM and ES. The left panels contain 
the mass and angular momentum accretion rates. The straight horizontal lines show the analytical mass accretion rates: dotted (panel a) is the 
Hoyle-Lyttleton rate (Eq. 1 in Paper I), solid are the Bondi-Hoyle approximation formula (Eq. 3 in Paper I) and half that value. The upper solid 
bold curve represents the numerically calculated mass accretion rate. The lower three curves of the left panels trace the x (solid), y (dashed) 
and z (dotted) component of the angular momentum accretion rate. The components are multiplied by the factor written down near the curves 
in order to amplify the deflections of the curves. The accretion rates of the components of the linear momentum (x solid, y dotted and z dashed) 
are plotted in the right panels. The components are multiplied by the factor written down near the curves in order to amplify the deflections of 
the curves. Additionally, the gravitational drag (x component) due to the asymmetrical mass distribution upstream and downstream from the 
accretor is shown as dash-dotted curve. 



lable 1. Parameters and some computed quantities tor all models. Moo is the Mach number of the unperturbed now, L size or the largest grid, 
R radius of the accretor, g = number of grid nesting depth levels, S = the size of one zone on the finest grid, e is the softening parameter (zones) 

for the potential of the accretor (see Paper I), U is the total time of the run, M is the integral average of the mass accretion rate, M is the 
maximum mass accretion rate, S is one standard deviation around the mean M of the mass accretion rate fluctuations, Mbh is denned in Eq. 3 
of Paper II, P is the approximate period of the accretion rate fluctuations ("vari": no periodicity is visible), t u is the approximate time at which 
the flow becomes unstable (the time at which the fluctuations of the ang. mom. accretion rate distincly rise over the background noise; this is 
the case when a value of 0.2 is reached), d ss shock standoff distance (measured from the center of the accretor; negative:downstream distance 
of the extrapolated shock lines), s is the entropy (Eq. 4 in Paper II); Z rms is the root mean square of the specific angular momentum magnitude; 
all time units are Ra/c^ ; all models have a specific heat ratio of 7 = 1 .01, the initial density distribution is randomly perturbed by 3%, and the 
number N of zones per grid dimension is 32. 
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3.1. Subsonic flow 

The distribution of matter for the three subsonic models can be 
seen in Fig. 1 , panels a, c and e. Since the accretor of model EL 
is large (one accretion radius), the gravity at and beyond the 
surface is weak. Thus upstream of the accretor matter gets 
absorbed faster than gravitational focusing can replenish the 
volume and so a lower density than at infinity is found around 
x = —2. Downstream of the accretor, at x = +2 a slight density 
enhancement is present due to gravitational focusing. For both 
models EM and ES, in which the accretor is much smaller (than 
model EL), gravitation can act more strongly and thus concen- 
tric density contours form. However, the center of the density 
contours does not coincide with the center of gravity which also 
is the center of the vacuum accreting volume. These features 
have been observed for the models with 7 = 4/3 and 7 = 5/3 
(cf. Fig. 1 in Papers III and IV). One very prominent difference 
between the models with different 7 is, however, the extent to 
which the accretor is eccentric relative to the concentric den- 
sity contours. While this shift is only marginal in the 7 = 5/3 
model SS (Paper III, Panel Id), it is approximately 0.5 R (R 
being the radius of the accretor) for the 7 = 4/3 model AS 
(Paper IV, Panel le), and for model ES (Fig. le) it is 1 R. The 
shift of the center of the contours is downstream relative to the 
accretor in all cases. 



Since matter is more compressible with decreasing values 
of 7, the maximum densities reached in the vicinity of the 
accretor increases when going from model SS over AS to ES. 
This behaviour can be seen when inspecting the size of some 
specific contour, e.g. lgp = 1, in the plots of these models. 

The accretion rates of several quantities for the subsonic 
models EL, EM and ES are shown in Fig. 2. A steady state is 
reached in all three cases with only small fluctuations in the 
rate of accreted linear and angular momentum. These are due 
to the initially perturbed medium and they die down after the 
inflowing unperturbed wind flows through the whole grid. A 
comparison of the accretion rates of these new models EL, EM, 
ES with models SL, SM, SS, AL, AM, AS from Papers III and 
IV reveal no striking differences. Only the sign of the accreted 
x-component of the linear momentum is clearly negative in 
model ES (Panel 2f) while it is clearly positive in model SS 
(Paper III, Fig. 2, bottom right). A positive sign indicates that 
more momentum is accreted from the forward (upwind) side of 
the accretor than from the back (downwind). 

Although it might seem that model ES is unstable, the fluc- 
tuations visible in Panels 2e and f are due to the initial condi- 
tions. It takes the flow roughly 27 time units (half the size of 
the grid divided by the flow speed) to cross the grid from the 
boundary to the accretor and thus to fill this area with 'fresh' 
unperturbed material. 



Since the models FL, GL and HL, simulating large accretors 
with size 1 .0 Ra, all show similar features, we will discuss 
these models together in this section. All these models reach a 
steady state in which the mass accretion rate is practically con- 
stant and the linear and angular momentum accretion rates are 
negligible, except, of course, for the x-component of the linear 
momentum (cf. Panels a and b in Figs. 3, 5 and 7). The value of 
the momentum's x-component is large and positive, indicating 
a large amount of momentum being accreted directly geomet- 
rically from the upstream flow. Model FL, however, does show 
fluctuating components of linear and angular momentum ac- 
cretion rates, albeit, at a very low level. The time it takes the 
initially perturbed density to cross half the grid is about 50 time 
units for model FL, so these fluctuations cannot be due to the 
initial perturbations. The density distribution and shock posi- 
tion at large distances (5 Ra or more) from the accretor are also 
representative for what the distribution looks like for the models 
with smaller accretor. They can be seen in Panels lb, and 4a 
and b. Obviously, as the Mach number increases the opening 
angle of the shock cone decreases (cf. Sect. 4.4). In all cases 
the shock is a static tail shock in which the flow downstream is 
directed away from the accretor. 

3.3. Mildly supersonic flow 

Contrary to the subsonic models, a shock front develops when 
the relative velocity between accretor and medium is super- 
sonic: contour plots of the density of the mildly supersonic 
models FL, FM and FS can be seen in Fig. 1 . In all these mod- 
els the accretor moves at a velocity of Mach 1.4 relative to 
the surrounding medium. Note the similarity of the contours 
of panels d and f in Fig. 1, manifest in the shape of the con- 
tours and the position and strength of the shock at some fixed 
distance. This indicates that the density distribution is not very 
sensitive to the size of the accretor for these models and to the 
differing mass accretion rates (see Table 1 and the following 
paragraphs) that they generate. It is due to the fact that the sonic 
distance is expected to be larger than the size of the accretor, as 
was mentioned in Section 2.2. Downstream of the accretor in 
models FM and FS (panels Id and f) a very spherically sym- 
metric flow develops within the shock cone. Contrary to the 
models BM, BS, KL and KS (in Papers III and IV, they have 
larger values of 7) the shock front in the models with 7 = 1 .01 
(panels Id and f) does not clear the surface of the accretor dur- 
ing the whole simulation. The fluctuations eventually die down 
and a steady state is reached. 

The large scale (around 10 Ra) density distribution and 
flow of model FL look very similar to those of the models with 
smaller accretor, while the flow at a distance of a few accretor 
radii from the surface of the accretor is somewhat different (see 
Fig. 1): Although the shock front in all three models ends at 
the surface and is a tail shock, practically no radial inflow is 
visible downstream for model FL (except within 0.3 accretor 
radii distance from the surface). 



the grid is about 10 time units for models FM and FS. Figure 3 
shows the runs of the accretion rates and of the gravitational 
drag as a function of time. The plot shows that these models 
converge toward a steady state after the initial transients have 
calmed down (the initial perturbations were advected off the 
grid). 

Also in the supersonic models FS and FM the x-component 
of the accreted linear momentum has a negative sign, indicating 
that more momentum is accreted from the back (downwind) side 
of the accretor than from the upwind side. Thus the accretion 
of linear momentum in these cases would not act as a drag but 
as an accelerating force! However, since the gravitational force 
always acts as drag and is stronger than the linear momentum 
accretion drag, the net force is decelerating. 

3.4. Moderately supersonic accretors 



density, model GS 
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Fig. 6. Gray scale plots showing the density distribution of model GS 
in a plane containing the center of the accretor. Darker shades represent 
higher densities. The time of the snapshot is the same as in Fig. 4e. 

Compared with all other simulations presented in this work, the 
moderately supersonic models GM and GS exhibit the most 
unstable flow and the largest fluctuations, e.g. the mass accre- 
tion rate varies by roughly 7% (cf. Table 1). This is, however, 
small compared to the models with larger 7 that produce mass 
accretion rate fluctuations of 14% (for 7 = 4/3) and of 23% 
(for 7 = 5/3). The comparison of different rates is done in 
more detail in Sect. 4. 1 . 

The contour plots of the density distribution (Fig. 4, pan- 
els c and e) show a striking difference to the equivalent plots 
of the models at different 7 (cf. Paper IV, Figs. 4c and 4e; Pa- 
per II, Figs. 6 and 11): the shock cone ends at the surface, i.e. it 
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Fig. 3. The accretion rates of several quantities are plotted as a function of time for all mildly supersonic models FL, FM and FS. The left panels 
contain the mass and angular momentum accretion rates, the right panels the linear momentum accretion rates and the gravitational drag. The 
components are multiplied by the factor written down near the curves in order to amplify the deflections of the curves. See also the caption of 
Fig. 2 
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Fig. 4. Contour plots showing snapshots of the density together with the flow pattern in a plane containing the center of the accretor for 
all moderately supersonic models GL (panel a), GM (panel c), GS (panel e) and highly supersonic models HL (panel b), HM (panel d), HS 
(panel f). The contour lines are spaced logarithmically in intervals of 0.1 dex for models GL, HL and GM, and 0.2 dex for the other models. 
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in the upper right hand corner of each three panels. 
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Fig. 5. The accretion rates of several quantities are plotted as a function of time for all moderately supersonic models GL, GM and GS. The 
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caption of Fig. 2 
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Fig. 7. The accretion rates of several quantities are plotted as a function of time for all highly supersonic models HL, HM and HS. The left 
panels contain the mass and angular momentum accretion rates, the right panels the linear momentum accretion rates and the gravitational drag. 
The components are multiplied by the factor written down near the curves in order to amplify the deflections of the curves. See also the caption 
of Fig. 2 



density distribution, one notices that the shock never clears the 
surface. This trend had already been apparent when comparing 
the results of the simulations for 7 = 5/3 to 7 = 4/3 (in Pa- 
per IV) - the shock front was closer to the accretor surface and 
touched the surface more often - and is now (with 7 = 1.01) 
even more extreme. So although the strength of the fluctuations 
is reduced in these models in which the shock ends at the ac- 
cretor surface, nevertheless, the flow is still unstable resulting 
in non-radial accretion flow. Lower density pockets, vortices 
etc. can be discerned. 

Two regions with a local density maximum are visible in 
Panel 4e and perhaps better visible in Fig. 6 : (a) directly down- 
stream in the wake of the accretor an accretion line forms and 
(b) directly along the shock cone after the infalling material has 
crossed the shock. It is similar to the mushroom shaped density 
distribution that has been described in Paper II, and that has 
preceded the unstable flow. In the new models presented in this 
paper the structure is fixed to the surface of the accretor, while 
for the 7 = 5 /3 flows the separate density maximum appeared 
when the shock was clearing the surface of the accretor. 

When inspecting the time evolution of the fluctuating ac- 
cretion rates of mass, angular and linear momentum (Fig. 5, 
Panels c-f) one sees the unstable flow reflected in the fluctuat- 
ing accretion rates. Obvious periodicities are not visible. 

3.5. Highly supersonic accretors 

The highly supersonic models HM and HS, too, do not con- 
verge to a quiescent steady state but show an unstable fluctuat- 
ing flow pattern. However, the fluctuations are less pronounced 
than in models GM and GS. In Panels 4d and f, snapshots of 
the typical density distribution have been plotted for the two 
models HM and HS, respectively. One cannot help but notice 
the very regular flow pattern and shock cone structure: the po- 
sition of the shock cone hardly changes (in time) and the shape 
is practically straight. Thus Panels d and f (of Fig. 4) just seem 
to be an enlargement of Panel b. Note, however, (a) the shock 
strength (more contour lines over a shorter distance) increases 
with decreasing size of the accretor and (b) there is a stagnation 
point at a position of about x = 0.8 (Panel 4d). Further out than 
this distance the flow escapes from the accretor (Panel b), while 
within this distance matter is accreted in a stream collimated 
toward the accretor (Panel f). This flow structure is very much 
different from the very complicated state shown in Fig. 8e in 
Paper IV (which is the equivalent model for 7 = 4/3). 

The temporal evolution of the values of some accreted quan- 
tities can be seen in Fig. 7. They show the slightly unstable 
accretion of mass and momenta. Contrary to the models GM 
and GS, the models HM and HS exhibit a periodicity with a 
period of approximately 0.2 time units, which can be seen from 
the Fourier transform of the mass accretion rates (Fig. 8). 

4. Analysis of Results 

4.1. Mass accretion rate 
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Fig. 8. The power spectrum vP{v) (in arbitrary units) of the mass 
accretion rate for models HM and HS. 
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Fig. 9. Mass accretion rates (units: 4nR^c c<D p ac ) as a function of ac- 
cretor size (units: Rb (Eq. 7 in Paper III), not accretion radii). The 
Mach number is an additional parameter and is specified by the num- 
bers to the right of the y-axis. The topmost bold curve corresponds 
to a Mach number of zero, i.e. a stationary accretor. The other lines 
belong to Mach numbers of 0.6, 1.4, 3.0 and 10., whereby the curve 
of Mach 0.6 completely overlaps with the curve belonging Mach 0. 
The numerical results to the corresponding Mach numbers are plotted 
with different symbols: the stars (*), triangles (A) crosses (x) and 
diamonds (o) are the results of models with A^oo=0.6, 1.4, 3.0 and 
10, respectively. The "error bars" extending from the symbols indicate 
two standard deviations from the mean (S in Table 1), for models in 
which the mass accretion rate fluctuates. 



Fig. 9. Note that in this figure the unit of length (for the accretor 
size) is .R B (Eq. 7 in Paper III), not Ra- In Paper III we had 
described an interpolation formula that had well fit (in the spirit 
of Bondi, 1952) the mass accretion rates found in Papers I, 
II and EH. Since all those models were done using a specific 
heat ratio of 5/3, the formula is adapted only to that one 7. In 
Paper IV the numeric mass accretion rates were compared to the 
values obtained with this formula and only a slight difference 
was found. In Fig. 9 we include the mass accretion rate that 
follows from our interpolation formula for 7 = 1.01 as solid 
curves. 

Two points seem noteworthy, (a) As has been mentioned 
in Sect. 2.2 the sonic distance for a gas with 7 = 1.01 is 
ci s = 0.49 R B . So just as was the case in Paper IV, the curves 
run horizontally for all accretors that have a radius smaller than 
the sonic distance. Within the sonic radius the flow is super- 
sonic and changes of the size of accretors smaller than the sonic 
distance should not influence the accretion flow. Within numer- 
ical accuracy, the mass accretion rates derived numerically are 
equal for those models in which the accretor size is smaller 
than the sonic distance, e.g. the mass accretion rate for mod- 
els HM and HS are 0.82 and 0.83, respectively (cf. Table 1). 
(b) The interpolation formula does not fit the mass accretion 
rates for all supersonic models too well: the numerical values 
lie above the interpolation curves by up to 60%. Especially the 
models with smaller Mach number deviate significantly from 
the formulae. 

The flow patterns are always less violent in the 7 = 1.01 
models compared to the 7 = 5/3 or 7 = 4/3 models: Fig. 10 
shows the relative mass fluctuations, i.e. thestandard deviation 
divided by the average mass accretion rate M/S (cf. Table 1), 
for all pertinent models. One notices that in all cases the largest 
bold symbols lie below the less large and bold symbols of 
same shape and same Mach number. This result extends the 
trend described in Paper IV refering to 7 = 4/3 compared to 
7 = 5/3 models. Again, the fluctuations of model HM are much 
smaller than those of model GM, just as had been the case with 
models DM and CM in Paper IV. It is due to the fact that in the 
models with high Mach number the shock hardly ever clears the 
surface of the accretor and when it does, it does not venture very 
far. Thus, the part of the accreted quantities that are accreted 
directly from the incoming stream (i.e. without having to cross 
the shock cone) remain constant and do not participate in the 
unstable accretion flow. 

The maximum density that is reached on the numerical grid 
is shown in Fig. 1 1 . In the same way as was described in Pa- 
per IV (and also in Shima et al., 1985) this density increases in 
models with ( 1 ) higher Mach numbers, (2) smaller accretors and 
(3) lower 7-values. Thus, these three points are only confirmed 
and extended in their trend by the new 7 = 1 .01 models. 

4.2. Angular momentum accretion rate 

For many applications the magnitude of the accreted specific 
angular momentum is of interest. In Fig. 12 we plot for several 
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Fig. 10. The relative mass fluctuations, i.e. the standard deviation S 
divided by the average mass accretion rate M (cf. Table 1), is shown 
as a function of the Mach number of the flow. Diamonds (o) denote 
models in which the accretor has a radius of 0.02 Ra, triangles (A) 
models with 0. 1 Ra- The largest bold symbols belong to models with 
7 = 1.01, while the smaller symbols belong to models with 7 = 4/3, 
and the smallest thin symbols are 7 = 5/3. The lines connecting the 
symbols are for reference only. 



models the ratio of temporal rms average of the specific angular 
momentum l rms over Mach number Mac as a function of the 
accretor size. The values of Z rms and can be found in 

Table 1. Additionally, at the top left of the plot, the diagonal 
line shows what values are expected analytically if one assumes 
that a vortex flows with Kepler velocity V around the accretor, 
just above the surface, i.e. Z rms is given by l s = RV = VRM X . 
This plot (Fig. 12) covers approximately the same range as 
Fig. 14 in Paper IV. Only those models that exhibit unstable, 
active flow patterns are shown. 

Note that the angular momentum accreted in the models 
presented in this paper (7 = 1.01) lie well below the values 
produced by the 7 = 5/3 models of Papers II and III, and by 
the 7 = 4/3 models of Paper IV. This confirms the fact that 
has already been mentioned: the fluctuations of the 7 = 1.01 
models are smaller than of the models with larger 7. 

4.3. Linear momentum accretion rate and drag forces 

Figures 13 and 14 are a collection of drag force values for all 
models that we have calculated upto now, from the series of 
Papers II to IV and the new results of the present paper. The 
amount of linear momentum accreted in all models is shown in 
Fig. 13. This linear momentum is what induces a drag force that 
can be called hydrodynamical drag. In the actual physical situ- 
ation the asymmetry of the mass density distribution upstream 
and downstream from the accretor produces an additional drag 




Fig. 11. The maximum density that is reached on the grid is shown as a 
function of the Mach number of the flow. Diamonds (o) denote models 
in which the accretor has a radius of 0.02 R A , triangles (A) models 
with 0.1 Ra, squares (□) models with I.ORa, while the cross (x) is 
model S6 with 0.2R A and the plus (+) is model T10 with 0.01i? A . 
The largest bold symbols belong to models with 7 = 1.01, while the 
smaller symbols belong to models with 7 = 4/3, and the smallest thin 
symbols are 7 = 5/3. The legend in the top left corner summarises 
how the nine symbols depend on the three values of 7 and on the three 
values of the accretor radius. The lines connecting the symbols are for 
reference only. 



force, that is called gravitational drag. This gravitational drag 
depends, however, logarithmically on the geometric extent of 
the surrounding medium through which the accretor moves 
(Chandrasekhar dynamical friction equations, e.g. in Binney 
& Tremaine, 1987; Chandrasekhar, 1943). For this reason the 
units of the gravitational drag values shown in Fig. 14 include 
the term ln(L / R A ), in which L is the size of the largest grid. 
One has to remember to remultiply this logarithmic term using 
the correct size of the medium when applying the plotted values 
to physical situations. 

A hydrodynamic drag force, contrary to an accelerating 
force, is actually only felt by the models simulating large accre- 
tors, with radii of one R A . In these cases, the geometric accretion 
of momentum from matter flowing directly from the upstream 
direction into the accretor (possibly slightly focussed by gravi- 
tation) dominates any other accretion from matter downstream. 
All subsonic models ( M^, = 0.6) produce hydrodynamic drag 
forces that are larger (more positive) than the corresponding 
supersonic models with other parameters kept equal. Thus the 
shock cone in the supersonic models seems to collimate the 
linear momentum accretion in the negative x-direction more 
strongly than the basically radial spherically symmetric but 
eccentric flow present in the subsonic models. Although the 
velocities are roughly of similar magnitude close to the accre- 



Fig. 12. As a function of accretor size the ratio of temporal rms average 
of the specific angular momentum Z rms over Mach number A4 00 is plot- 
ted for all models: squares (□) for models FS and FM from Paper III, 
crosses ( x ) for the models from Paper II, diamonds (o) for models CM 
and CS from Paper IV, triangles (A) for models DM and DS from Pa- 
per IV, large bold stars (*) for models GM and GS and finally large 
bold plus-signs (+) for models HM and HS. The diagonal straight line 
at the top left shows what values Z s (see text) are obtained if the flow 
would rotate with Kepler velocity along the accretor surface. 

tor, independent of whether the matter has passed the shock or 
not, the density is larger within the shock cone, so the momen- 
tum accreted from the shock cone is larger than the momentum 
accreted from the upstream direction. For this reason the mo- 
mentum accreted by the supersonic models has a larger share of 
negative momentum. Thus most models of interest (supersonic 
and not large) actually experience an accelerating momentum 
accretion force. Only the smallest subsonic models (ES and EM) 
have a slightly negative accreted momentum, too, generated by 
the very eccentric density contours (cf. Fig. 1), that produce a 
high density just downstream of the accretor. The gravitational 
force always acts as drag (decelerating) and is stronger than the 
momentum accretion drag in all our models, so the net force 
is decelerating. However, it seems in principle possible to con- 
struct a situation such that the gravitational drag is smaller than 
the magnitude of the hydrodynamic drag. The latter can be neg- 
ative, so a net acceleration would result. Of course, energy and 
momentum conservation dictates, that this situation cannot be 
a stationary one. 

Another feature visible in Fig. 13 is that with decreasing 
7 the accreted momentum becomes more negative for models 
with accretor sizes 0.02 R A and 0.1 R A . This, too, can be traced 
to the more regular and better collimated flow in the models 
with smaller 7, which can directly be seen when comparing the 
contour plots, e.g. Fig. 4e and f with Fig. 6e and Fig. 8e from 
Paper IV. 
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Fig. 13. The linear momentum accretion rate is plotted as a function of 
accretor size for all models : the 7 = 5/3 values are taken from Papers II 
and III, the 7 = 4/3 values from Paper IV and the 7 = 1.01 values are 
showninFigs. 2,3, 5, and 7. The bold curve over the horizontal straight 
line (separating positive and negative values) represents the amount of 
momentum that would be accreted purely through the geometric cross 
section (i.e. -k ^pv^). The units are Po = nRxpv 2 ^; negative values 
indicate an acceleration, positive values a drag. 



Fig. 14. The gravitational drag force density is plotted as a function of 
accretor size for all models: the 7 = 5/3 values are taken from Papers II 
and III, the 7 = 4/3 values from Paper IV and the 7 = 1.01 values are 
shown in Figs. 2, 3, 5, and 7. The units are Fo = Po Iii(L/Ra), with 
Po defined in the caption of Fig. 13, and L is the size of the largest 
grid (cf. Table 1). The zero line is plotted for reference only. 



or negligible (within numeric accuracy) gravitational drag. Just 
as with the hydrodynamic drag (Fig. 13) the spread of values 
due to different Mach numbers of the flow is of the same or- 
der of magnitude as the spread due to different accretor sizes. 
A clear trend can be seen, that (a) with decreasing accretor 
sizes (b) for smaller Mach numbers and (c) for smaller 7, the 
gravitational drag force increses. The softer equation of state 
in models with smaller values of 7 yields larger accretion cone 
densities and thus stronger gravitational forces. Although flows 
with higher Mach numbers produce stronger shocks and thus 
higher densities within the shock cone, too, the shock openening 
angle in these fast flows is smaller. Thus the gain in gravita- 
tional drag due to the higher density seems to be offset by the 
smaller shock cone. The density increases close to the accretor, 
so it is not surprising that models with smaller accretor sizes 
show a larger gravitational drag: the front-back asymmetry is 
larger. Note the big difference in magnitude of the gravitational 
drag when comparing the subsonic models (Mach 0.6) with the 
mildly supersonic models (Mach 1.4): the drag increases by 
an order of magnitude. A further increase of the Mach number 
from 1.6 to 10, on the other hand, reduces by roughly 30% the 
gravitational drag (when assuming a fixed accretor size). Thus 
the presence of the shock cone is crucial for a non-negligible 
gravitational drag force. 



4.4. Comparison of shock opening angles 
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Fig. 15. The shock opening angles for various models are plotted versus 
speed of the flow. The vertical lines show the numerical models, solid 
represents the 7 = 1.01 models, dashed the 7=4/3 models and 
dotted the 7 = 5/3 models. The vertical lines have slightly been 
shifted horizontally to facilitate the distinction of the lines. The curve 
decreasing shows the relation = arcsin(l/.M). 



per IV), the shock opening angle (8, angle of the shock front 
to the direction of the unperturbed flow) seen in the simula- 
tions tends to be larger than the values predicted analytically, 
9 = arcsin(l/A4), for large distances from the accretor. A 
rough estimate of the shock opening angle can be obtained 
from the contour plots and the values are plotted in Fig. 15 by 
vertical bars, the range of which span the different angles of the 
inner and outer shock cone edges, and of the different models, 
at large distances from the accretor. 

Contrary to the results presented in Paper IV for 7 = 4/3, 
the deviation of the numerically obtained angles from the ana- 
lytic curve for models with 7 = 1 .01 is about equal for all Mach 
numbers. The suggestion by other authors (e.g. Livio et al., 
1991) that it is the pressure within the shock cone that enlarges 
and expands the boundaries of the cone (thus increasing the 
opening angle) seems to be corroborated: for the models with 
a large Mach number (Moo = 10) the shock angle decreases 
with decreasing 7. Smaller 7 imply a softer equation of state 
which potentially produces a smaller pressure. That this effect 
is strongest at large Mach numbers is due to the strong shock 
and thus high density and pressure generated in the wake. 

5. Conclusions 

We draw the following conclusions from the results presented 
here: 

1. For 7 = 1.01 the flow (and thus the accretion rate of all 
variables) exhibits active unstable phases only for models 
in which the flow speed is large enough and the accretor 
size small enough: all models with Mach numbers of 0.6 
or 1.4 tend toward a steady state with quiescent accretion 
flow, as well as the Mach 3 and 10 models with an accretor 
size of li?A- Only the models at Mach 3 and 10 and with 
accretor sizes 0. 1 and 0.02 show unstable flow patterns. 

2. A periodicity with period 0.2 time units can be discerned 
for the Mach 10 models. 

3. The shock front is always attached to the accretor and forms 
a tail shock, contrary to most simulations of flow with 7 = 
5 /3 and 7 = 4/3. The fluctuations generated in the unstable 
models are smaller in the isothermal simulations than in the 
models with larger 7. 

4. The smaller value of 7 (softer equation of state) has as 
consequence, that higher densities have to be reached before 
a sufficiently high pressure is built up to support the shock 
front (compared to models with larger 7). Thus the densities 
in the shock cone are higher and still the pressure is not 
enough to separate the shock from the accretor (tail shock). 
Both these consequences then produce a higher rate on 
linear momentum accretion and a smaller fluctuation of the 
cone and accretion rates. 

5. The model with Mach 0.6 yields accretion rates identical (to 
numerical accuracy) to those in stationary accretor models. 

6. The mass accretion rates are slightly larger for the 7 = 1 .01 
models compared to 7 = 4/3. The phenomenological fit we 
proposed in a previous paper and which worked well for the 



obtained values by as much as 60% for 7 = 1 .01 . A&A style file L-AA version 3. 

7. The rms of the specific angular momentum accreted in the 
isothermal simulations is about a factor of 2-3 smaller than 
in the models for 7 = 5 /3 and 7 = 4/3. 

8. A comparison of the hydrodynamical drag (= accreted linear 
momentum) yields: (1) only the accretors with radius 1 R A 
are subjected to significant drag. Smaller accretors in sub- 
sonic models experience hardly any drag at all, while small 
accretors in supersonic models are accelerated (i.e. the rel- 
ative velocity between accretor and medium is increased) 
by the momentum they accrete. However, the gravitational 
force always acts as drag and is stronger than the momen- 
tum accretion drag, so the net force is decelerating. (2) The 
acceleration is larger with decreasing 7. 

9. A comparison of the gravitational drag (front-back 
density asymmetry) yields: (1) it is always positive, 
i.e. decelerating. (2) The subsonic models experience a drag 
roughly one order of magnitude smaller than the supersonic 
models. (3) The drag decreases with increasing accretor 
size, with increasing Mach number and with increasing 7. 

Movies in mpeg format of the dynamical evolution 
of all models are available in the WWW at http : 
/ / www . mpa-garching . mpg . de/~ mor/bhla . html 
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